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We report that upon excitation by a single pulse, the classical harmonic oscillator immersed in 
classical electromagnetic zero-point radiation, as described by random electrodynamics, exhibits 
a quantized excitation spectrum in agreement to that of the quantum harmonic oscillator. This 
numerical result is interesting in view of the generally accepted idea that classical theories do not 
support quantized energy spectra. 
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Historically, the discreteness of atomic spectra moti- 
vated the early development of quantum mechanics. Be- 
cause quantum mechanics has made numerous success- 
ful predictions on atomic spectra, conventionally the dis- 
crete atomic spectrum is attributed to transitions be- 
tween the quantized energy levels of an atom. Classi- 
cal theories, such as stochastic electrodynamics (SED), 
are usually thought of as incapable of generating discrete 
atomic spectra, or discrete spectra for any potential that 
supports quantum mechanical bound states. Indeed, Pe- 
ter Milonni has commented in his well-known book The 
Quantum Vacuum that "Being a purely classical theory 

of radiation and matter, SED is unable to account 

for the discrete energy levels of the interacting atoms." 

IB- 

In the case of the harmonic oscillator, quantum me- 
chanics predicts a discrete multi-resonance spectrum at 
the harmonics of the oscillator's natural frequency. The 
classical harmonic oscillator only supports a single reso- 
nance at the oscillator's natural frequency. In this Let- 
ter, we investigate the excitation spectrum of a one- 
dimensional harmonic oscillator immersed in the classi- 
cal zero-point field, or zero-point radiation. The classical 
zero-point field (ZPF) is defined through the theory of 
random electrodynamic^] [2j [3 , as a nonzero homogenous 
solution to Maxwell's equations. The Plank constant ft is 
brought into the classical theory as an overall factor that 
sets the strength of the classical zero-point field. Our sim- 
ulation shows that a classical harmonic oscillator in the 
classical ZPF exhibits the same discrete multi-resonance 
excitation spectrum as a quantum harmonic oscillator. 

The excitation spectra of the classical and the quantum 
system are obtained through numerical computation of 
the total energy of the oscillator after excitation by a sin- 
gle pulse. The quantum mechanical result is obtained by 
numerically solving Schrodinger's equation with a high- 
order multipole expansion. 
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1 Random electrodynamics is also known as stochastic electrody- 
namics. 



First, we discuss the theory for a one-dimensional 
quantum harmonic oscillator. The Hamiltonian for an 
oscillator interacting with a pulse is 
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where = (^,0,0), x = (£,0,0), q is the charge, m 
is the mass, ujo is the natural frequency, and ((/> P ,A P ) 
is the potential of the pulse. Using the Coulomb gauge 
V • A p = 0, Maxwell's equations yield qb p = in the 
absence of external charge. In addition, p x -A p = A^-p^. 
Therefore, the Hamiltonian becomes 
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The eigenvalues of the first part of the Hamiltonian, 
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give the discrete energy levels E n = huJo(n-\-l/2) = huj n 
for the unperturbed quantum harmonic oscillator. The 
oscillator is initially in its ground state. The second part 
of the Hamiltonian, 



H =-^-(2A p -p x -qA 2 p ), 
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can induce transitions to the higher excited states. The 
field A p is described by a propagating Gaussian pulse 



A p = Aq cos (k p • x — Upt) exp 
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where At is the temporal width of the pulse, Ax = cAt is 
the spatial width, k p = u p /c (sin pi 0, cos 6 P ) is the wave 
vector of the carrier wave, and s p = (cos# p , 0, — sm0 p ) 
is the polarization vector of the pulse. The field ampli- 
tude is Aq. After excitation, the state \^)(t) becomes a 
superposition of the energy eigenstates |n), 
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FIG. 1. The sampled vacuum modes in k-space. Left: The sampled vacuum modes (black dots) are distributed in k-space as 
a spherical shell with thickness A/c. The number of sampled modes shown is N u — 2 x 10 4 . A slice of the spherical shell at 
k x ~ is highlighted (red dots). Right: A quarter of the highlighted slice (red dots) is projected on the k y k z -plane. The modes 
sampled in cartesian coordinates (red circles) are shown for comparison. 



To obtain the coefficients, c n , of the excited state 
|^)(t), we solve the Schrodinger equation, 
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where C(t) is a TV x 1 matrix and the Hamiltonian H 
is represented by a N x TV matrix, 



the complex field 
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is approximated by 
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For each matrix element in the matrix of H , uj nrn = 
u n —ujm and H nrn = (n\H \m). To approximate a system 
with an infinite number of energy levels N we choose 
N = 20 for numerical convergence. Since resonances at 
the harmonics of the oscillator's natural frequency are 
affected by the spatial dependence of the field, the dipole 
approximation k • x ~ 1 is not sufficient for the study of 
the excitation spectrum, making a numerical approach to 
this problem convenient. In our simulation, we expand 
the spatial dependence of the field A p = (^A p + Aj^ /2 
up to the 20th order. Thus, the spatial dependence of 
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Lastly, the matrix element of the operators, x and p x , 
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FIG. 2. The excitation spectra of harmonic oscillators in different theories. The energy absorption spectrum is shown for the 
harmonic oscillator described in quantum mechanics (left top), classical ZPF theory (left middle), and classical mechanics (left 
bottom). The results from quantum mechanics and classical ZPF theory (i.e random electrodynamics) are in agreement. The 
vertical axis shows the ensemble averaged energy for the two classical theories and the expectation value of the energy for 
quantum mechanics. The energy is scaled by Tiujq. The horizontal axis shows the pulse frequency uj v scaled by the oscillator's 
natural frequency ujq. The natural frequency is chosen to be c^o = 10 16 (rad/s), the charge is chosen to be q = 1.60 x 10 _19 (C), 
and the mass is chosen to be m = 9.11 x 10 _35 (kg). The choice of the mass is made to keep the integration time manageable 
without losing the physical characteristics of the problem. The ZPF frequency bandwidth is set as A = 2.2 x 10 2 Tcj 2 ,, which 
is much larger than the resonance width of a classical oscillator (i.e. Tuq). The calculations are convergent for increasing ZPF 
frequency bandwidth. The number of sampled ZPF modes is N u = 500. Although this number of modes is much smaller than 
the one shown in FIG.^ N u = 500 is sufficient for reaching numerical convergence. 



are specified by 
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where e x = (1, 0, 0). 

Next, we discuss the theory for an one-dimensional 
classical harmonic oscillator in the classical ZPF [3l H] . 
The classical Hamiltonian corresponding to the quantum 
system, Eq. is 
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The classical ZPF in the Hamiltonian, H c i, is specified 
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where uj = c|k|, # k A is the random phase uniformly 
distributed in [0,27r], and V is the volume of a cavity. 
The two unit vectors, e k and s k 2 , describe a mutually 
orthogonal polarization basis in a plane perpendicular to 
the wave vector k. The pulse field A p is identical to that 
in Eq. ([5| except for x being a classical quantity rather 
than an operator, 



A p = Aq cos (k p • x — ujpt) exp 
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From the Hamiltonian, H c i, the classical equation of 
motion can be derived, 
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FIG. 3. Excitation spectra of a harmonic oscillator at different excitation pulse angles P . The excitation spectra of a harmonic 
oscillator in classical ZPF theory (red solid line) are compared to those in quantum mechanics (blue dash line). For different 
excitation pulse angles P = 0, 7r/6, 7r/4, tt/3, 7r/2, the spectra is shown. When the pulse is perpendicular to the direction of 
the oscillator's motion (i.e. 9 P = 0), there is no higher harmonic resonance. When the pulse is parallel to the direction of the 
oscillator's motion (i.e. 6 P — 7r/2), there is no excitation at all. 



where mT'x is the phenomenological radiation damping 
term added to the theory. The damping coefficient is 
^_ 2q 2 1 
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and the magnetic field by 



B = V x A. The symbol denotes the ^-component 
of the vector E. Because v x = (^,0,0), the magnetic 
part of the Lorentz force is zero. To avoid numerical 
runaway solutions, we assume the point-particle descrip- 
tion of the charged particle and make the approximation 
mT'x ~ — mTuj^x [6j[7]. Hence, the equation of motion 
becomes 
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which will be used in our numerical simulation. 

To carry out the simulation, a set of modes (k^, A) 
is chosen to construct the classical ZPF in Eq. (16) [4]. 
The wave vectors Iq have frequencies within the finite 
range [oj - A/2, oj + A/2], where A is the ZPF fre- 
quency bandwidth and is much larger than the oscilla- 
tor's resonance width TcUq. The ZPF modes in k-space 
are sampled in spherical coordinates. Such a sampling 
method is computationally more efficient in reaching nu- 
merical convergence and approaches sampling in carte- 
sian coordinates as — » oo (see FIG.[T]). The specific 
steps of the sampling method are given in the following. 
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and 



Ki = (uj - A/2) 3 /3c 3 + (i - 1)Ak 

* = 

The stepsize An is specified by 

[(w„ + A/2) 3 /3c 3 - (w - A/2) 3 /3c 3 ] 



(22) 



Ak ■ 



N,„ 



(23) 



The random number is uniformly distributed in 

[—1,1], and is uniformly distributed in [0, 2tt]. Fi- 
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nally, the polarization vectors are evaluated according to 
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For a large number of modes, TV^, the sampling in spher- 
ical coordinates and cartesian coordinates approach each 



other, and the volume factor V in Eq. (16) can be esti- 
mated by 



atomic 1/r potential alone affords the possibility of sub- 
harmonic excitation [9], in our work the nonlinearity is 
due to the existence of the classical ZPF and the spatial 
dependence of the excitation pulse, as exemplified by the 
necessity of calculation beyond the dipole approximation. 
Our numerical approach may also be useful to critically 
test recent claims that superposition and entanglement 
are supported by classical ZPF theories [lOj [11] . 
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Here, the k-space volume encloses the sampled ZPF 
modes Iq in a spherical shell, 




The result of our numerical simulation for the excita- 
tion spectrum is given in FIG. [2] When the classical ZPF 
is absent, the classical harmonic oscillator supports only 
a singe resonance at its natural frequency. With the clas- 
sical ZPF acting as a constant background perturbation, 
the classical harmonic oscillator exhibits a discrete multi- 
resonance excitation spectrum with the background en- 
ergy shifted up to Tiuoq/2 [3, 4 . The shape and the mag- 
nitude of the resonance peaks are in agreement with the 
quantum mechanical result. Furthermore, as the pulse 
angle 6 P changes from to 7r/2, the energy of the har- 
monic oscillator in the classical ZPF scales in the same 
way as that of the quantum harmonic oscillator, as shown 
in FIG. H 

In conclusion, we have shown that the classical har- 
monic oscillator in the classical ZPF displays the same 
discrete multi-resonance excitation spectrum as the 
quantum harmonic oscillator. The agreement between 
classical ZPF theory and quantum mechanics is within 
the error of our numerical convergence. Extension of nu- 
merical work to atomic system appears to be interest- 
ing [8]. Whereas in Cole's work the nonlinearity in the 
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